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ABSTRACT 



We investigate a model QCD sum rule for the pion wave function Lp-„{x) based 
on the non-diagonal correlator whose perturbative spectral density vanishes and 
, ^{x, M^), the theoretical side of the sum rule, consists of condensate contributions 

' only. We study the dependence of <i>(x, M"^) on the Borel parameter and observe 

, that $(a;,M^) has a humpy form, with the humps becoming more and more pro- 

\^ ' nounced when Af ^ increases. We demonstrate that this phenomenon reflects just 

, the oscillatory nature of the higher states wave functions, while the lowest state 

' wave function extracted from our QCD sum rule analysis, has no humps, is 

rather narrow and its shape is close to the asymptotic form ip'^'^{x) = 6x(l — x). 



CIh' 1. QCD sum rules and pion wave function 



The pion wave function (p-rrix) is the basic object in the perturbative QCD (pQCD) description 
of hard exclusive processes involving the pion: (j)Trix) is the probability amplitude to find the pion 
' in a state composed of its two valence quarks carrying the fractions xP and (1 — x)P of its large 

■ longitudinal momentum P. More rigorously, the pion wave function (pT^{x,fi) can be defined as the 

function, whose A^-th moment is given by the matrix element of a local operator with N covariant 
derivatives Bi: 

{P^'P"^ ...P"''} [ (t){x;^x)x^ dx = i^~^{Q\d-i5{-i''D''^ ...D''''}u\ti+,P)\^ (1) 

"'0 

where {. . .} denotes the symmetric-traceless part of a tensor and /i is the renormalization parameter 
for the composite operator On- Instead of x, it is more convenient sometimes to use the relative 
variable ^ defined by a; = (1 + ^)/2. In the limit of exact u-d symmetry, the pion wave function is 
an even function of ^, i.e., all odd ^-moments of 4>{x') A*) vanish. This definition of the wave function 
implies that its integral normalization is fixed by the matrix element of the axial current: 

= /tt, (2) 
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Figure 1: Normalized pion wave function ipT^{x): asymptotic limit (solid 
line) and CZ-model (dashed line). 



where f^r ~ 133 MeV is the pion decay constant. In many cases, it is convenient to use the normalized 
wave function </?^(a;;^) = 0^(a;;/i)//^ whose integral is simply 1. The /^-dependence of lPt^{x; ^) is 
governed by the evolution equation which follows from the renormalization-group equation for the 
composite operator J^. In the /.t — > oo limit (p-n-{x, ^) has a simple and natural form i (see alsoi) 



However, one is usually interested in the form of (p^rix; fj,) at low values of the renormalization 
parameter /i ~ 1 GeV relevant to experimentally accessible situations. This form is determined 
by non-pcrturbative QCD dynamics and, in principle, it may strongly differ from the asymptotic 
limit. To calculate the pion wave function at low values of the probing parameter /i, one should take 
into account non-perturbative aspects of QCD. The closest to pQCD and one of the most popular 
non-perturbative approaches is provided by QCD sum rule£ which incorporate information about 
the non-trivial structure of the QCD vacuum via the operator product expansion. 

The first application of QCD sum rules to the pion wave function (p-nix) was the calculation of 
the pion decay constant i.e., the zeroth moment of 4>-k{x), in the pioneering paper by Shifman, 
Vainshtein and Zakharov □ who considered the correlator of two axial currents and calculated 
within 5% accuracy. Next step was made by Chernyak and A.Zhitnitskyi who tried to construct 
the whole pion wave function lyS^r {x) using information about its lowest moments calculated from the 
QCD sum rule analysis for correlators of the axial current with the above mentioned local operators 
On. The model wave function chosen by Czi 



has the moments {i"^)^^ ~ 0.43 and {S,'^)^^ « 0.28, to be compared with {^'^Y" = 0.2 and = 
3/35 for the asymptotic wave function. The large value of these moments dictates the characteristic 
double-humped shape of the CZ wave function: it has maxima at a: ~ 0.15 and x f» 0.85 and a 
zero at the middlepoint x = 0.5 (see Fig.l). Such a form for the wave function of a lowest state 
looks rather strange. Quantum-mechanics-based intuition would rather suggest that the ground 
state wave function has a shape like that of the asymptotic wave function: without humps, nodes or 
zeros. One can also expect that all such peculiarities should appear - and in increasing number(!) 
- for the wave functions of radial excitations. 



fTr{x] /i OO) = <p$^'*(x) — 6x(l — x). 



(3) 



if'^^{x) = 30x(l - x){\ - 2xf, 



(4) 



Formally, the reason for the exotic shape of the CZ model wave function can be traced to the 
structure of the nonperturbative (condensate) terms in their sum rule. Being written directly for 
the wave function (^^(x), it reads: 

f^ifi^ix) + {higher states) = - x) + "l^f^l [Six) +6(1- x)] 

+ ^^^^^^{imx) + 6il-x)]+2[S'ix)+S'il-x)]}. (5) 

The perturbative loop contribution in this sum rule has a smooth behavior coinciding with the 
asymptotic a;(l — x) shape. On the other hand, the condensate terms are strongly peaked at the 
end-points x — and x — 1, i.e., in the regions where one of the quarks has zero momentum. The 
standard ansatz for the higher states is to model them by perturbative spectral density, which is 
~ a;(l — x) in this case. As a result, there remains only one state, the pion, whose wave function has 
to reflect the presence of the condensate peaks at a; = and x = 1. In other words, the CZ wave 
function looks like a compromise between the smooth perturbative loop behavior ~ x{l — x) and 
the condensate contributions forcing strong enhancements at a; = and x — I. 

Earheii, we argued that taking the CZ sum rule (^) at face value amounts to assumption that 
vacuum quarks have zero momentum, which is an approximation with a limited applicability range. 
In general, one would expect that vacuum quarks have a smooth distribution in momentum, and, 
hence, a d{x) term, say, should be treated only as the first term of an expansion of a smooth function 
in a series over S(x) and its derivatives. If the generating smooth function is not very narrow, then 
the condensate peaks are not as drastic as the S{x) + 6(1 — x)-approximation, and the impact of the 
condensate corrections on the pion wave function is much milderQ. 



2. Nondiagonal correlator 



In what follows, we would like to concentrate on another subtle point of the standard QCD 
sum rule analysis of hadronic wave functions, namely, on the implicit assumption that higher states 
can be modeled by ("are dual to") a perturbative spectral density. In fact, this assumption is in 
an obvious conflict with a standard quantum-mechanical situation, when the ground state has a 
monotonous positive-definite wave function, its first radial excitation has one zero, the second has 
two and higher state wave functions become more and more oscillating. To study this problem in 
its cleanest form, it makes sense to analize a sum rule with vanishing perturbative density. This can 
be easily arranged by taking a non-diagonal correlatoiH, e.g., the correlator of the generic operator 
On with the pseudoscalar current dj^u rather than with the axial current d'y^^nU = Oq. Then, for 
massless quarks, all the perturbative terms (i.e., those corresponding to the unity operator) of the 
operator product expansion are zero because an odd number of gamma-matrices would be involved 
in any trace, and only condensate terms appear on the theoretical side of the sum rule. To the lowest 
order of a^, all these terms have a singular behaviour like 6{x), S'(x), etc. 

The correlator of d^^u and u'y^^D^ d has a remarkably simple structure in the limiting case 
= both on the theoretical and phenomenological sides of the sum rule: 

a) there is only one particle ~ the pion - which has nonzero projections on both the axial current 
d'y5Ji_iU and the pseudoscalar current d'y^u ; 

b) for massless quarks, a single operator - the quark condensate - survives in the operator product 
expansion (a proof can be found in the SVZ paper§). Two other, formally allowed terms (qD^q) 



and {q{aG)q) cancel each other (recall that (qD^q) ~ ^{qig{aG)q)) . This leads to the well-known 
PCAC relation 

(0|J75w|7r) = ^{{uu) + {dd)). (6) 

The relevant QCD sum rule, in its borelized form, with being the Borel parameter charac- 
terizing the exponential suppression of the higher states contribution, looks as follows: 

i^). + {e).'e-<-"'' + {e)."e-<"'^'" + [Mgher states) = i±i_ll!^ + 0{l/M% (7) 

where 0{\/NP) includes the power suppressed contributions due to higher condensates {qD^q), 
{q{(jG)q), etc. These contributions vanish for = 0, i.e., their coefficients contain factor N. 

Furthermore, since only the pion term survives on the l.h.s. of this sum rule in the specific case 
= 0, we have: 

= (C^=")., = (e"'=")..=... = 0. (8) 

A simple observation is that the ("axial") wave functions of the higher pseudoscalar mesons tt', 
tt", . . . must (!) have oscillations to produce zero total integrals. In principle, the pion wave function 
may have oscillations or humps as well, but this is not mandatory. 
As before, we rewrite the sum rule directly for the wave functions: 

^.(x) + ^.,ix)e-<'/''' +...= ^('')+^^(^-'^) + a{qD'q) {s' (x) + 5'(1 ^ x)} + . . . . (9) 

Note, that higher condensates cannot have 6{x) or (5(1 — x) coefficients, since all higher condensate 
terms must disappear after one takes the zeroth moment of this sum rule. 

Our last comment here is that the smooth functions ipjr{x),ipTr'{x), . . . on the l.h.s. of the 
sum rule can be produced only by an infinite summation of singular distributions (5"(x),(5"(l — x) 
associated with the local condensates. 

3. Nonlocal condensates 

In the coordinate representation, the contribution of the simplest diagram is given by the product 
of the perturbative propagator S{z) ~ (z^7^)/2;^ and the nonlocal condensate {q{0)q{z)). Next 
term (evaluated in the Fock-Schwinger gauge, which is the most convenient for the QCD sum rule 
calculations) is proportional to {q{0){(jG{0))q{z)) {z'^j^)/z'^. Performing the Taylor expansion of the 
nonlocal condensates in produces the OPE in terms of local condensates. The resulting sum rule 
for the wave function has the structure of an expansion over the delta functions ^(a:), (5(1 — x) and 
their derivatives. 

Our strategy0 is to avoid the Taylor expansion to preserve the smoothness properties of the 
objects involved on the theoretical side of the sum rule, i.e., keep together all terms generated by a 
particular nonlocal condensate. As the next step, we construct model expressions for the nonlocal 
condensates to see how the properties of the nonlocal condensates affect the form of the pion wave 
function extracted from the relevant sum rule. 

It is convenient to parametrize the z^-dependence of the simplest bilocal quark condensate 
{q{0)q{z)) = {q{0)q{0))Q{z^) with the help of a Laplace-type representation 

/"OO 

Q(z2)= / e'^'/^f{s)ds. (10) 
Jo 



The spectral function /(s) may be called "the distribution function of quarks in the vacuum" since 
its Tith moment is proportional to the matrix element of the local operator with to rtth power: 







In particular, for the lowest two moments one has 



(11) 



and 







fis)ds = l (12) 



with having the meaning of the average virtuality of vacuum quarks. 

In a similar way, parametrizing the quark-gluon nonlocal condensate {q(0)ig{crG{0))q{z)) = 
{qig{aG)q)Qi{z^), one can introduce the quark-gluon distribution function fi{s). Since 

{qD\) = ]^{qig{<jG)q), (14) 
there exist a relation between the zeroth moment of fi{s) and the first moment of /(s): 

ml= / /i(s)ds = 4 / sf{s)ds. (15) 
Jo Jo 

The standard QCD sum rule estimatci for ttiq is TOq ~ 0.8 GeF^; since — rtv^j^^ one can take 
\\ = OAGeV^. 

Constructing models of nonlocal condensates, one should satisfy also some other constraints. 
For instance, if one assumes that the vacuum matrix element {q{D'^)^° q) exists, then /(s) should 
vanish faster than l/s^°+^ as s oo. So, if all such matrix elements exist, /(s) must vanish faster 
than any power of 1/s for large s. As a possible choice, one may impose that, at large s, the function 
f{s) behaves like /(s) ~ e^**^/""^ (Gaussian fall-off), or /(s) ~ e^'^l'^ (exponential fall-off), etc. The 
opposite, small-s limit of /(s) is governed by the large-|z| properties of the function Q(z^), i.e., by 
its behaviour at large space separations or at large values of the imaginary time variable r — izQ. 
The latter case can be easily assessed using the QCD sum rule for the heavy-light meson spectrum 
in the heavy quark effective theory (HQET). In the HQET, the heavy quark has a trivial propagator 
Sq{z) ^ d^{z)d{zQ) and, hence, the time dependence of the correlator of two heavy-light currents is 
determined by the light quark propagator!. At large imaginary time r, the correlator is dominated 
by the lowest state contribution ^ e""^^ where A = (Mg — mQ)\mQ^<xi is the lowest energy level of 
the mesons in HQET. This means that Q{z^) ^ e"'^'"^ for large Euclidean z and /(s) ~ e~^^/'* in 
the small-s region. Numerically, A is around 0.45 GeV. 

Combining, in the simplest way, the e^^^^" dependence with, say, the Gaussian fall-off at large 
s, we arrive at the ansatz 

f{s) = Ne-'f^V'^-sya^ ^^g) 

where A^ = 0.2 GeV^, the normalization constant N is fixed by eq.(p^) and the cr-parameter is fixed 
by eq.(|l^, where for the average virtuality of vacuum quarks we take the usual QCD sum rule 
value! Xl ~ OAGeV^. 





{GeV'-2) 



Figure 2: Our model for the vacuum distribution function /(s). 

4. Sum rule 

Now, the sum rule for the pion wave function can be written as 

(^.(x) + ^^,{x)e-<'l''" +^^.{x)e-<"'^'" +... = ^x,AP) (17) 
— (1 - x)f{xKP) + -h(xM^) + F{x, M2)| + {x^l-x) 

where /, /i , . . . are the vacuum distribution functions corresponding to the lowest nonlocal con- 
densates {qq) and {qaGq), respectively, and the function F{x,M^) is given by higher nonlocal 
condensates {qGGq), (qGGGq) ,etc. 

As discussed above, the sum rule reduces to an extremely simple form if one takes its zeroth 
moment: 

1 

<^{x,AP)dx = 1, (18) 



with only the simplest quark condensate term contributing to "1" in the r.h.s. This imposes the 
following relations for the two lowest terms of the nonlocal condensate expansion: 

f{s) ds - 1, 





h{s)ds = A sf{s)ds. (19) 
10 Jo 

The second relation reflects the cancellation between {qaGq) and (qD^q) terms of the local expansion. 
For remaining condensates we have 



/ F{x,]\'P)dx ^0. 
Jo 



(20) 



By virtue of eqs.([T^, ^5|), the first two relations are satisfied for = oo and since, at large s, the 
distribution functions are supposed to vanish faster than any power of 1/s, the violation of the above 
finite-M^ relations also drops faster than any power of 1/M^ at large M^. This is consistent with 
the fact that contributions decreasing faster than any power may be missed by the operator product 



expansion. The practical lesson is that the Borel parameter should be taken in the region where 
the violation of the normalization condition ( p^ ) is sufficiently small. 

To fix the sum rule, one should specify a model for the quark-gluon nonlocal condensate and 
say something about the higher contributions denoted by F{x,M'^). As we will see, our procedure 
of extracting wave functions from the sum rule is perfectly linear, in the sense that each nonlocal 
condensate term on the theoretical side of the sum rule produces an additive contribution to all the 
wave functions on its l.h.s. Hence, one can split each wave function into respective parts generated 
by a) the lowest quark condensate, 6) quark-gluon condensate, c) next condensate, etc. Then one 
can study separately the resulting sum rules, each having only one type of the nonlocal condensate 
on its r.h.s. Finally, one should add the contributions extracted from each of these partial sum rules. 
In fact, to illustrate general features of the fitting procedure, it is sufficient to analize the sum rule 
containing the lowest nonlocal condensate on its theoretical side. However, since the two lowest 
nonlocal condensates are related by eq. ([l9|), these two terms should be better considered together. 

The structure of the nonlocal quark-gluon condensate is specified by the relevant distribution 
function /i(s). Information about this function, in principle, can be also obtained from a (future) 
study of the QCD sum rules in the heavy quark limit. Lacking such information at the moment, 
we will assume the simplified ansatz that /i(s) coincides with the function /(s) governing the 2^- 
dependence of the simplest nonlocal quark condensate. This assumption is not crucial and it does 
not affect qualitative features of our analysis. 



5. Fitting sum rule 



Now we can write down our model sum rule for the "axial" wave functions of the pseudoscalar 
mesons: 

^^{x) + ^^,{x)e-<'/^'" +^^,,{x)e~<"^^'" + ... = ^x,M^) (21) 

= J,)/(.An+(.^i-x) 

with the function /(s) specified in the preceding section. For the Tr'-meson, we will take the experi- 
mental mass TO J, ~ 1.7 GeV^. 

It is evident from this sum rule that the function $(a;,M^), i.e., the weighted sum of all wave 
functions is given by two humps, which are moving as changes. As Af^ increases, the humps 
become narrower, higher and move towards respective boundary points x ^ or x = I, approaching 
the S{x) or d{l — x) form in the Af^ 00 limit. For AP — 1 GeV^, e.g., the function ^{x, M^) looks 
very much like the Chernyak-Zhitnitsky wave function (see Fig. 3). However, one should remember 
that ^{x, M^) is not just equal to the pion wave function: at large there might be a large 
contamination from higher states. Taking larger Af^, e.g., Af^ ~ 1.2GeV^ produces even a wider 
function, while decreasing Af^ to 0.8 GeV^^ produces a function with closer and lower humps. 

The lower Af^, the more pronounced is the dominance of the pion in the total sum ^{x,M'^). 
However, we cannot take too low M^, because the operator product expansion might fail. Since the 
average virtuality of the vacuum quarks is O.iGeV'^, it is definitely unreasonable to go below 
the point M'^ = O.AGeV^, because our "large" probing virtuality M'^ should be larger than - 
otherwise one should expand the correlator in 1/A^ rather than in 1/Af^. Taking A4^ — 0.4GeV^, 



Figure 3: Function M^) for = 1 GeV^ (left), = 1.2 GeV^ 
(middle) and = O.SGeV^ (right). 



we observe that $(a;, M^) is very close to the asymptotic wave function of the pion (see Fig. 4). 
Assuming that the total sum $(a;, M^) at such low is completely dominated by the pion, we 
have to conclude that our model for the nonlocal condensate sum rule suggests that the pion wave 
function is rather close to its asymptotic form. However, one should be more accurate here, since 
even the modest increases of to 0.5 GeV^ or 0.6 GeV^ induce humps in ^{x,M^) (see Fig. 4). 




Figure 4: Fmiction $(a;, M^) for M ^ = 0.4 GeV^ (left), = 0.5 GeV"^ 
(middle) and = 0.6 GeV^ (right). 

This means that the Tr'-contribution is visible at '--^ 0.5 GeV"^, and one should better try to 
fit <I>(x, M^) by two lowest states. Taking two different but close values of M^, one can extract the 
relevant wave functions, 

^{x,Ml)e<'^^'" - $(.T,M|)e<'/^2 

Choosing — 0.5 Ge^^ and M\ = 0.55 GeV^ we obtained the curve for ^Pt^{x) shown in Fig. 5 
(left). A very close result is obtained if one takes the pair Mf = 0.55 GeV^ and = 0.60 GeF^ 
(see Fig.5, middle). For definiteness, we will fix the pion wave function as that extracted from the 
first pair Ml = 0.5 GeV^ and = 0.55 GeF^. The relevant n' wave function is then shown in 
Fig.5 (right). Note, that its maxima are by a factor of 10 higher than those of fix). However, 
this only means that the overall scale characterizing the magnitude of matrix elements (0| . . . |7r') is 
essentially larger than But this is only natural in view of large mass of the tt' particle. 

To estimate the contribution of the resonances higher than n', it is convenient to introduce the 



Figure 5: Pion wave function extracted from the two-states fit of 
$(x, M2) performed a) at = 0.5 GeV"^ and M| = 0.55 GeV^ (left) 
and b) at Mf = 0.55 GeV^ and M| = 0.6 GeV^ (middle), c) Wave 
function of the 7r'-mcson extracted from the two-states fit of $(a;, M^) 
performed at = 0.5 GeV^ and M| = 0.55 GeV^ (right). 

function 

xAx,M^) = mx,M^) - ip^ix)) e<'l^\ (23) 

At low M^, this function is very close to the tt' wave function ip-K'ix) as determined from our two- 
states fit. In particular, for the reference points = 0.5 GeV^ and = 0.55 GeF^, this function 
simply coincides with (/)^'(.t). For larger Af^, however, the higher resonances modify its form more 
and more strongly. This can be seen from Fig. 6. It is evident that the difference between x-n\^^ M"^^ 
and ^-Kiix) = XTr'{x,M'^ = 0.55 GeV^) increases with M^. One can guess that the increase of 
X%'{x,M^) — ip„'{x) just reflects the increasing contribution of the next resonance. Looking at the 
actual curves for the difference x-k' {x, M^) - cp^, (x) at three values, = 0.8 GeV'^ , A'P = 1 GeV"^ 
and ~ 1.2 GeV^, (see Fig. 6 (right)), one can notice that, to good accuracy, the difference 
XTr'{x,M^) — (Pt,'{x) has essentially the same shape for different M^, with absolute normalization 
governed by an M^-dependent factor. 




Figure 6: Function x^'{x,M'^) for = 0.8 GeV^ (left), and = 
IGeV^ (middle), shown together with (p.jr'{x) (dashed line). Right pic- 
ture shows the increase of the difference XTr'ix,]^"^) — (fv'{x) when the 
Borel parameter takes the values = 0.8, 1 and 1.2 GeV^. 

Now, the question is whether one can fit the combination 

^x,M^) - M^) - ^A^)e-<'l^" = {Xn'{x,M^) - ip^,{x))e-<'^^' (24) 



Figure 7: Mass difference mj^ 
x-regions. 



i'^, calculated via eq.(26) in essential 



by the next resonance contribution ipii{x)e '"bA^ , This means that we should take the ratio 



Xu'(x,M|) 
and try to see whether it can be fitted by 



(25) 



^,nl,/M? -ml/Ml 
m?^,IMl-myMl 



This task can be reformulated as a procedure determining the mass of the third resonance from the 
relation 



rriR 



Again, we take two pairs: a) {M^, Af|) 



In 



(1, 0.8) GeV^ and b) {Mf, M|) 



(26) 

(1.2, 0.8) GeV^ and 

plot the r.h.s. of eq.(p6|) for three x-regions (see Fig. 7). Because of the zeros of Xtv' (x, M^) — lPt^i {x), 
the curves are not as constant as one might wish. Still, one can safely state that m|j — m^, = 
2.5 ± 0.5 GeV^ is a reasonable estimate. For the mass itself, this gives rn|j = 4.2 ± 0.5 GeV^ ^ which 
translates into a rather narrow prediction for the tt" mass: m^" = 2.05 ± 0.15 GeV. 




Figure 8: Pion wave functions from Fig. 4 (left) and their difference 
(right). 

Though the third resonance is rather massive, it is not completely invisible in the low-M^ region. 
In particular, looking at the difference between the two pion wave functions shown in Fig. 4 (recall 
that they were extracted from two-states fits performed for different pairs) , one can see that the 
resulting curve has the shape specific for the third resonance (Fig. 8). 



6. Summary and conclusions 

In this paper, we considered a model sum rule for the pion wave function. A specific feature of 
this sum rule is that the usual perturbative contribution is absent altogether, and the theoretical 
side of the sum rule is given by the condensate contributions only. To represent the latter, we 
incorporated nonlocal condensates. As a result, $(.t,M^), the weighted sum of wave fimctions 
related to pion and its radial excitations, was given by a curve generated by two humps which were 
moving to the end points x = and x = 1, raising in height with increasing Borel parameter M^. On 
the other hand, the relative weight of the higher states increases when gets larger. This clearly 
indicates that the peaks observed in <I>(x, M^) for > 0.4 GeV"^ reflect only the oscillatory nature 
of the wave functions related to the pion excitations. Our explicit fits confirmed this expectation: 
the pion wave function extracted from this sum rule has no humps and is rather narrow, despite all 
the humpy nature of $(a;, M^). This result has evident implications for the CZ sum rule based on a 
diagonal correlator: the peaks in the end-point regions generated by the condensates contribution, 
reflect the oscillatory components in the higher states wave functions rather than the humpy wave 
function of the ground state, the pion. 
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